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In this paper, we study power allocation strategies for a fixed-gain amplify-and-forward relay network employing 
multiple relays. We consider two optimization problems for the relay network: 1) optimal power allocation to 



, maximize the end-to-end signal-to-noise ratio (SNR) and 2) minimizing the total consumed power while maintaining 

> 

the end-to-end SNR over a threshold value. We investigate these two problems for two relaying protocols of ali- 
en 1 ' 

On : participate relaying and selective relaying and multiple cases of available channel state information (CSI) at the 

relays. We show that the SNR maximization problem is concave and the power minimization problem is convex for 
all protocols and CSI cases considered. We obtain closed-form expressions for the two problems in the case for full 
CSI and CSI of all the relay-destination links at the relays and solve the problems through convex programming 
when full CSI or CSI of the relay-destination links are not available at the relays. Numerical results show the 
benefit of having full CSI at the relays for both optimization problems. However, they also show that CSI overhead 
can be reduced by having only partial CSI at the relays with only a small degradation in performance. 
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I. Introduction 

The classical three-terminal relay channel has been around since the 1960s when it was first introduced 
by Van Der Meulen [1], [2]. Cover and Gamal then characterized the capacity of a three-terminal 
relay channel where a relay assists a source node in communicating with a destination node [3]. Since 
then precious little work was carried out on relays until a decade ago. The rapid progress in wireless 
communications technology, the increasing popularity of tetherless connectivity and satisfying the high 
quality-of-service (QoS) requirements have rekindled interest in relays [4], [5]. Recent years have seen a 
dearth of work being carried out to study the performance of relay-assisted systems [6]— [13]. It has been 
shown that relays can provide spatial diversity [14]— [16], increase capacity [17]— [19], conserve power 
[20], [21] and enhance coverage [22]-[24]. 

There are two main protocols for relays: 1) amplify-and-forward (AF) (non-regenerative) relays, 2) 
decode-and-forward (DF) (regenerative relays) [25]. AF relays amplify the received signal and forward it 
to the intended destination. AF relays are unable to detect any errors in the received signal and hence, 
are termed as non-regenerative. Unlike AF relays, DF relays can detect the presence of an error and 
can correct it before re-transmission and hence, are termed as regenerative. However, DF relays require 
extra processing that increases complexity at the relays which often needs to be simple. In the following, 
we consider AF relays which can be further divided into two groups; namely, fixed-gain AF relays and 
variable-gain AF relays [26]. In fixed-gain AF relays, the relay gain does not depend on the fading gain 
of the source-relay link and remains constant throughout the duration of system operation, whereas in 
variable-gain AF relays, the relay gain changes with the fading gain of the source-relay link. Thus, in 
variable-gain relays, the relay always requires CSI of the source-relay link. In this work, we concentrate 
on fixed-gain AF relays. 

Power is a precious resource. It was reported in [27] that information and communication technology 
(ICT) utilizes more than 3% of the total electrical energy consumed worldwide and this percentage is 
expected to increase with time. Furthermore, in wireless communications, vendors and operators are 
searching for energy-efficient algorithms and devices to cut down on energy and operating costs [28]. 
In addition, mobile devices need to conserve energy as they have limited battery. Therefore, it is of 
paramount importance to use the available power as efficiently as possible. In a cooperative relay network, 
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this corresponds to allocating power efficiently among the source and relay nodes to improve performance. 
Moreover, it is also essential to reduce overhead in the system. Most of the works available on power 
allocation for relay networks focus on the case of full CSI available at the relays. This usually requires 
the controller 1 to feedback CSI of all the links to each relay. This causes great overhead and consumes 
precious system resources such as bandwidth and time. Thus, it is crucial to come up with power allocation 
strategies which can work with partial CSI or channel statistics at the relay to minimize the feedback 
overhead [29]. However, this reduced overhead comes at the cost of decrease in performance as there is 
less knowledge to work with and exploit. Therefore, there is a performance- overhead trade-off. It is good 
to have insight into this trade-off, so that optimal decisions can be taken to improve system performance 
under different scenarios. 

Power allocation for cooperative network with fixed-gain AF relays was studied in [30]-[34] 2 . In [30], 
the authors proposed power allocation schemes to maximize the sum and product of the average SNR 
of the source-destination and average SNR of the relay-destination link for a fixed gain relay-assisted 
source-destination pair. Moreover, the schemes required only the knowledge of channel statistics. For the 
same model as in [30], [31] proposed optimal and near-optimal power allocation algorithms to maximize 
the end-to-end SNR under slowly varying channel conditions. In [32], power allocation to minimize the 
symbol error rate (SER) for M-ary phase shift keying (MPSK) was derived for a source communicating 
with destination through a single fixed-gain AF relay. It was shown in [32] that the power allocation 
method provided better performance when the relay was near the destination. References [33] and [34], 
both proposed power allocation strategies to minimize the bit error rate (BER) of a communication system 
employing a single fixed-gain AF relay. In [35], we studied power allocations for an AP fixed gain AF 
relay network with complete CSI of the relay-destination link and knowledge of only channel statistics 
of the source-destination and source-relay links. We provided closed-form solutions for the optimization 
problems considered 3 . 

In this work, we consider a source-destination pair which communicates with the help of m fixed-gain 
AF relays. In addition to the m dual-hop links, the source and destination are also connected through a 

'The place where system decisions are taken, for instance it can be the destination. 

2 There are many other works on power allocation for relays, however most of them focus on variable-gain relays. 

3 We made a slight mistake for the implementation of the energy-efficiency problem solution in [35] which makes the results presented in 
[35] a little worse than they actually should be. This mistake and correction will be discussed in Section II-C1 



4 

direct link 4 . The relays are assumed to have only one antenna and work in half-duplex mode. Thus, they 
cannot transmit and receive at the same time and at the same frequency. To avoid interference, all the 
relays are assumed to operate on orthogonal channels 5 . For this system, we consider two relay participation 
schemes of all-participate (AP) relaying in which all the m relays forward the signal to the destination 
and selective relaying in which only the selected 6 relay forwards the signal received from the source to 
the destination. For both schemes, we consider two optimization problems of optimal power allocation 
(OPA) and energy-efficiency. We refer to these as dual problems. OPA refers to the problem of allocating 
power to the source and the relays to maximize the end-to-end SNR under a total power constraint on the 
system. In the dual problem of energy-efficiency, the total power consumed is minimized while keeping 
the end-to-end SNR above a certain threshold. Furthermore, we consider three cases of available CSI at 
the relays. In all the three cases, it is assumed that the destination has full CSI. The three cases of CSI 
considered are: 

1) Full CSI of all the links at the relays. 

2) Knowledge of the channel statistics of all the links at the relays. 

3) Two cases of partial CSI at the relays 

a) Full CSI of all the relay-destination links and knowledge of only channel statistics of all the 
source-relay and source-destination links (Partial CSI-/3). 

b) Full CSI of all the source-relay and source-destination links and knowledge of only channel 
statistics of all the relay-destination links (Partial CSI-a). 

We study the two dual optimization problems for the two relay participation schemes under all three cases 
of CSI given above. 

Our main contribution in this paper is to give a detailed analysis of the system under AP relaying, 
selective relaying and the three CSI cases. We give insight on the performance-overhead trade-off. We show 
which system configurations perform better for which objectives under which conditions. For instance, 
numerical results show that when the fading gains of the source-relay links and the source-destination 
links are the same on average, then having knowledge of the instantaneous CSI of the source-relay links 

4 The results in this paper include the case of no direct link as a special scenario by replacing the fading gain of the source-destination 
link by 0. 

5 They can be orthogonal in time, frequency, space, or code. 

6 How selection takes place is discussed in Section III when we consider selective relaying in detail. 
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benefits more than having knowledge of the source-destination links for OPA. 

The remainder of the paper is organized as follows. In Section II, we consider AP relaying. We discuss 
both optimization problems of OPA and energy-efficiency for all the CSI cases listed above. Section III 
focuses on selective relaying. We again discuss OPA and energy-efficiency under all three assumptions on 
the CSI. Numerical results are presented in Section IV. Finally, Section V summarizes the main results 
of the paper. 

II. All-Participate Scheme 

A. System Model 

Consider a cooperative system with a source node, a destination node, and m relays. Each relay is 
assumed to be equipped with only one antenna and operates in half-duplex mode. The source and the 
relays transmit on orthogonal channels. Without loss of generality, we assume time division multiple access 
(TDMA). The transmission takes place in two phases. In the first phase, the source transmits information 
to the relays and the destination. The signals at the ith relay and the destination are given by 

y si = \fE~h si s + n si (1) 

y sd = y/E~h sd s + n sd , (2) 

where E s is the source energy, h si and h sd are the instantaneous channel gains from the source to the ith 
relay and destination respectively. The channels are modeled as independent Nakagami-m random vari- 
ables. n si and n id are the complex Gaussian noise with zero mean and variances <j 2 si and a 2 sd , respectively. 
In the second phase, the relays, after amplification, forward the signal to the destination. The received 
signal at the destination from the ith relay is 

y id = \fa~E~E i h si h id s + n { , (3) 

where a { is the ith relay gain, E,- L is the ith relay energy, h id is the instantaneous channel gain from the 
ith relay to the destination which is again modeled as Nakagami-m fading, n,i is rij ~ CiV(0, of) and 

2 7-1 I t 1 2 2 2 

a i = aihji\h id \ a si + a id . 
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One can write the m + 1 received signals at the destination in matrix form 

y = hs + n, (4) 



where 



y = 



11 l 
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a 2 sd y OiEi | /iid 1 2 < + <r 2 d y a lJ E m | /imd I V s 2 m + < d 

and n ~ CN (0,1). Throughout this paper, it is assumed that the destination has complete CSI of all the 
links. It is also assumed that all the links experience independent fading. Furthermore, the fading gain of 
each link changes independently from one time slot to another. 

B. Optimal Power Allocation 

First, we consider the problem of OPA. In OPA, the end-to-end SNR is maximized under power 
constraints on the system. For OPA, we consider the three cases of full CSI at the relays, knowledge of 
only channel statistics at the relays and partial CSI at the relays. 

I) Full CSI: In this section, we assume that the relays also have complete CSI of all the links. Using 
(4), and assuming maximal-ratio-combining (MRC) at the destination, the end-to-end SNR of the system 
is given by 

/ m m j. \ 



,i=0 



where a = Mr-, a* = , = -4r- and Q — \. In this work, we consider a total power constraint 

^sd °"si ® id ° si 



on the whole system and individual power constraints on all the nodes. Hence, the power allocation 
problem can be written as 

max 7, subject to 

Es,Ei 1 

(6) 

E s + J2 E i< E tot, 0<E S < E™\ 0<Ei< E™ ax . 

In general, 7 is not a concave function of the source and relay powers as its Hessian is not negative 
semi-definite. However, as we show in Appendix A, the objective function in (6) is concave given the 



constraints on the system. Hence, the optimization problem, (6), can be solved using the Lagrange dual 
method [36]. Moreover, as the constraints are affine, Slater's condition [36] is satisfied, i.e. the duality 
gap between the primal and dual solutions is zero. Therefore, the solution obtained for the Lagrange dual 
problem is also the optimal solution of the primal problem in (6). Hence, solving the problem in (6) using 
the Lagrange dual method, with the help of [35], yields the solution 
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From equations (7) and (8), one can conclude that the OPA follows a water-filling solution. Hence, power 
is allocated in an iterative manner. However, unlike traditional water-filling algorithm, here the closed- 
form solution may change according to the initial results. If the source power, E s , satisfies its constraints 
and any relay power does not satisfy its individual constraint, then the solution to the problem remains the 
same, however, the optimization variables and the constraint changes. Therefore, in this case, the optimal 
solution is given by 
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where X, Y and Z represent the sets of powers which satisfy the individual constraints, are less than zero 
and greater than the peak individual constraint, respectively. Now, if the source power is greater than 
E™ ax , then it is set at E™ ax . The updated optimal solution for the relay powers now becomes 

rpmax 
Hi ^ 



where the Lagrangian multiplier can be obtained as 

rpmax I / a j Cj \ 



o 



Et0t ^ + ^^J 3 ) 

iez jex / 
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2) Knowledge of only Channel Statistics: In Section II-B1, it was assumed that the each relay had 
full CSI of all the links. However, this greatly increases the complexity of the system. The destination 
has to inform all relays of the CSI through dedicated feedback channels. This consumes a considerable 
amount of resources. Moreover, the reverse link between the destination and the relays might be poor 
and communication might not be possible. Therefore, it is desirable to be able to work with less CSI. 
Hence, in this section, we assume that each relay only has knowledge of the channel statistics of all the 
links and the destination has full CSI. Thus, to perform power allocation, the end-to-end SNR needs to 
be averaged over all the links. As the channels are modeled as Nakagami-m random variable, ctj and $ 
are both Gamma random variables and their probability density functions and are given by 

fa-(x) = — t— x^-V^ and f .(x) = l -— x k ^~ 1 e~^l x > 0, (13) 

rftj^ r(fc A )7jf 

respectively, where k ai and are the shape parameters of the links, % z and 7^ are the average SNRs of 
the links and is the gamma function [37, Eq. (8.310.1)] As all the links are assumed to be independent, 



the average end-to-end SNR can be found by averaging (5) over the density functions in (13) 



(14) 
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Solving the above with the help of [37, Eq. (3.383.10)] gives the average end-to-end SNR as 

7 = E. £ - E, £ *^ e =rt*r (i - k 0l , -JL-) (15) 

where .) is the upper incomplete gamma function given by T(s, z) = J 2 °° e'H^dt [37, Eq. (8.350.2)]. 
The optimization problem is the same as in II-B1, however the objective function is now changed to 7 
instead of 7. We show that 7 is a concave function of the optimization parameters on the domain of 
interest in Appendix B. Thus, the optimization problem is concave. However, it is difficult to find closed- 
form expressions for the optimal solution due to the complexity of the objective function. Fortunately, 
as the problem is concave, we can utilize well-known algorithms for convex optimization. So, one such 
algorithm, the interior point algorithm can be used to find the optimal solution [36]. 
For the special case of Rayleigh fading, 7 in (15) simplifies to 

where Ei(.) is the exponential integral function of the first kind and is related to the exponential integral 
function as E ± (x) = -E^-x) [38]. 

3) Partial CSI: Two important sub-cases of the full CSI and knowledge of channel statistics cases 
discussed in the previous two sections are the partial CSI instances: 1) partial CSI-/3 2) partial CSI-a. 
These are important as they give us an idea regarding the trade-off between performance and complexity. 
Studying the performance of these two cases of partial CSI provides us insight as to knowledge of which 
link is more important. Thus, if an increase in performance is desired, then only the links which have 
greater affect on the performance of the system can be fed back to the relays. 

For the first case, the objective function can be obtained from 7 by replacing a, L by k a .% r Then power 
allocation is done according to the algorithm described in Section II-B1. In the second case, the objective 
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function can be obtained from 7 by replacing A; a .7 a . by a { . Then resource allocation can be performed 
using the interior-point algorithm as in Section II-B2. 

C. Energy-Efficiency 

In Section II-B, we considered the problem of maximizing the end-to-end SNR under individual and 
total power constraints. In this section, we study the problem of energy-efficiency. The objective is to 
minimize the total power consumed, E tot = E s + Y1T=\ ^i, while keeping the instantaneous end-to-end 
SNR above a threshold, r f h , and ensuring the source and relay powers do not exceed their respective 
individual constraints. Like Section II-B, we consider the three cases of full CSI, knowledge of only 
channel statistics and partial CSI. 

1 ) Full CSI: With full CSI at the relay, the optimization problem is given by 



Now, we need to show that the optimization problem in (17) is convex and hence, the optimal solution 
can be found. The objective function and the individual constraints are convex function. So, it is only 
left to show that 7 is concave, meaning ^ th — 7 is convex, on the domain of the problem. Employing 
the same notation and procedure as in Appendix A, to prove concavity of 7, we need to show that 
D(x,y) = (f(x) - f(y))(g(x) - g(y)) < 0. If f(x) > f(y), then g(y) > g(x), to satisfy the constraint 
on 7 and vice versa. Thus, D(x,y) < and 7 is concave on the domain. Moreover, 7 is a monotonically 
increasing function of E s and E^s, hence the optimal solution to (17) is achieved when 7 = ^ th . As 
the other two constraints are affine and the objective function is convex, Slater's condition is satisfied. 
Therefore, the solution obtained using the Lagrange dual method will be optimal. Solving the problem 
(17) using the Lagrange dual method, with the help of [35], gives the optimal solution as 



min E t 

Es,Ei 



subject to 



(17) 



7 > 7 , < E s < E™ ax , 0<Ei< El 



max 
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(18) 








(19) 
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where 



P 



\^ =1 V W 



+ 



(20) 



From (18) and (19), one can see that the solution again follows a water-filling algorithm described in 
Section II-B1. Hence, the power allocation process is repeated until all the powers satisfy the constraint. 
However, the optimal solution changes depending on the initial power allocation. There are two cases like 
for the problem of OPA: 1) E s lies between ans E™ ax , 2) E s is greater than E r s nax . Considering case 1 
first, the optimal power allocation is given by 
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(21) 



As stated above, the power allocation in this section follows the same procedure as in Section II-B1. 
Hence, the power is allocated in an iterative manner. However, the channel conditions can be such that 
the solution obtained is where some powers are set at their maximum constraints and some at 0. A 
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situation can arise in this instance where 7 > j th due to the powers greater than their constraint being 
treated first in the algorithm. Hence, there should be a check at the end of the algorithm and if 7 > j th , 
then the whole water-filling process is repeated, however, this time the powers which came out to be zero 
in the first run of the complete procedure are not included in the algorithm from the start 7 . 

2) Knowledge of only Channel Statistics: Now, we consider the case of knowledge of only statistics 
at the relays. As was the case in Section II-B2, the end-to-end SNR now needs to be averaged over all 
the links. Therefore, in this instance, the energy-efficiency optimization problem is given by 

min E tot , subject to 

Es,Ei (22) 

7 > 7 t/l , < E s < E™ ax , 0<Ei< E™ ax , 

where 7 is given by (15). Using Appendix B and a similar argument as that in Section II-C1, it can 
be shown that the the optimization problem, (22), is convex. Hence, the Lagrange multiplier and other 
convex optimization algorithms can be applied to solve the problem. However, due to the complexity of 
the problem, it is very difficult to obtain closed-form expressions for the optimal solution. Therefore, we 
use the interior-point algorithm as utilized in Section II-B2 to obtain the optimal solution. 

3) Partial CSI: In this section, we study the the energy-efficiency problem for the two cases of partial 
CSI. As was the case in Section II-B3, the both the partial CSI scenarios are special cases of the full CSI 
case. The objective functions and the power allocation algorithms can be obtained from the Section II-C1 
for partial CSI-/3 and Section II-C2 for partial CSI-a as was the case in Section II-B3 

III. Selection Scheme 

In Section II, AP relaying was considered. However, AP relaying requires additional complexity at the 
destination to combine the relays. Also, as the relays transmit on orthogonal channels, it consumes a huge 
amount of system resources and decreases throughput. To ameliorate these drawbacks of AP relaying, 
selective relaying has been proposed in which only the "best" relay is selected to forward the signal from 
the source to the destination. The selection criteria depends on the objective. For OPA, the relay which 
maximizes the end-to-end SNR after power allocation is selected. For energy-efficiency, the relay which 
minimizes the consumed energy while fulfilling the constraint on the end-to-end SNR is selected. If no 

7 This was our mistake in [35]. We did not account for these special cases 
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relay fulfills the constraint on the end-to-end SNR, then the relay which achieves the maximum end-to-end 
SNR is selected. 



A. Optimal Power Allocation 

1 ) Full CSI: For selective relaying, noting that the power is now only divided between the source and 
one relay, one can re-formulate the end-to-end SNR in (5), in the case that the ith relay is selected to 
transmit, as 

7i = i]i E tot (oq + a, t a : 7 ) i = 1, 2, . . . , m, (23) 

where we have replaced E s = rjiE tot , Ei — (1 — r\i)E tot and < r\i < 1. Ignoring the individual power 
constraints, the optimization problem is 

max riiEtot ( a + an ^ - - ) i = 1, 2, . . . , to. (24) 

The concavity of the objective function follows from the concavity of the problem for AP relaying. Taking 
the derivative of the objective function in (24) and equating it to zero yields the optimal solution 



n 1 / / (diEtotPi + QaiCi . 

where and ^ are the links associated with the ith relay. If ^ is found to be such that one of the powers 
exceeds its individual constraints, then r\i is adjusted that the power lies on its peak individual constraint. 
The case where both powers exceed their constraints is when E tot > E™ ax + E™ ax . In this case, both the 
source and the selected relay transmit at their individual constraints. The algorithm for power allocation 
for selective relaying is 

• Calculate rji for all the relays using (25). 

• Compute the resulting % for each relay. 

• Select the relay which has the maximum 7,. 

2) Knowledge of only Channel Statistics: As was the case for AP relaying, if CSI is not available at the 
relays, then the end-to-end SNR has to be averaged over the channels before power allocation is performed. 
Hence, now the relay which gives best performance on average is now selected. The optimization problem 
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now, with the help of (15), is 

max (1 - m )E tot L 0%0 + k ao%i - ka %£ 1 e^F (l - k(k , — JM ) (26) 

where now E s = (1 — r)i)E to t, Ei = rjiE to t and < < 1 due to ease of analysis. The concavity of 
the problem follows directly from the concavity of the AP case. Taking the derivative of the objective 
function in (26) and equating it to zero gives 

a^^-^V ^^ +2 T " ft 'a^^J e (27) 



*ah^m r ( 1 - fc Aj — A— ) + 



Equation (27) can be solved numerically using algorithms such as bisection, Newton's method etc to yield 
the optimal value of r/j. Similar to the full CSI case, r)i is found for all the relays and then the relay which 
maximizes the averaged end-to-end SNR is selected. For the special case of Rayleigh fading, rji can be 
obtained from 

= ~E tot U ao + - lai - e^^lE, ( ^— ) ) - (1 - r/)% 

_}_ e ^Am i Ei ( ] + —(—-E 1 ( ^— ) e awE%tot ^- 



Vf \aiViEtotlfsJ Vi \Vi \aiViEtotlpJ a^E^p. 

3) Partial CSI: The two partial CSI cases can again be handled in the same way as in Section II-B3. 
The objective functions and the power allocation algorithms can be obtained from the full CSI case in 
Section III-A1 and the channel statistics case in Section III-A2 in the manner as outlined in Section II-B3. 

B. Energy-Efficiency 

I ) Full CSI: In the case of full CSI, the energy-efficiency problem for the ith selected relay is 

min E tot , subject to 

Es,Ei (29) 

7i > 1 th , 0<E S < E™ ax , 0<Ei< E™ ax , 
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7j . Bs|Q0 + ai _^L, (30 ) 



where 

— K_ I ren 4- rv„- — 

OiEipi + Q 

The optimization problem, (29), is solved for all the relays and the relay which minimizes E tot while 
fulfilling the constraint on ^ is selected. If no relay fulfills the constraint on 7$, then the relay which 
maximizes 7$ is selected. 

Ignoring the individual constraints and taking advantage of the fact that at the optimal solution 7 = 7^, 
we can write 

aiEiBa th + Ca th 

E s = — — — . (31) 

a Q + aiEiaoPi + diEiaiPi 

Thus, ignoring the individual constraints, (29) can be re-written as 

mm h iij. (32) 

E i aoQ + aiEiaoPi + diEiCtifii 

Taking the derivative and equating it to gives as 



Substituting (33) in (31) yields 



hi = in — ; — ^ — • ' 



e, = w^+Mt,. (34) 

(ce + oii) a/ aiOti/3i(a th 



Incorporating the individual constraints gives the water- filling solution 

gmax 



, , . (Vaj(XiPiti + C»Qq) 1 th . MVs 



(a + di) V CLiCtiftiCa 



th 







y/ aiaif3i(a th - aoCi \ ' q 6) 
ciiPi (a + ai ) J 

2) Knowledge of only Channel Statistics: In this case, the selection procedure and the optimization 
problem are the same as in Section III-B1, however the constraint on the end-to-end SNR changes to 
7i > 1 th , where 

h 'v C k<h 1 c -i f t- 



***** V Pl, a t E af3i 
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Again exploiting the equality of the on j t , we obtain 



th E k Pi 

E s = \ B c , . (38) 

(k ao % + vyoj - **^ e -*^* r (1 - 

Therefore, the optimization problem becomes 

min 1 \ g c . + ^. (39) 

(Wy«„ + k a a ai ) E^ - k -^4^e^ r (l - 

Taking the derivative and equating to gives 

^ fc a;7o:;Sj „ ,-B,*7fl. r> / 1 7. Si \ I 7 ^Oi^ftT^Sj 



(40) 

The above equation can be solved through bisection search to yield the value of Ei which can be substituted 
back into (38) to obtain E s . The maximum of E s and Ei is checked and if it exceeds its peak constraint, 
then it is set at it its peak constraint and the other power is obtained from the constraint. If no power 
exceeds its respective peak constraint, then the minimum power is checked and if it is below 0, it is set 
to zero and the other power is obtained from the constraint. 
For Rayleigh fading, (41) simplifies to 

(7a + Taj Ei - ^ e 7 ^E 1 ( —^—) ) ' - f^^e^l E 1 ^ C 



th^ A th^ A f a \ u 

a iEfli3 t diErfp. l \a i Ea p J 
3) Partial CSI: The two partial CSI cases can be obtained from the full CSI case and the channel 
statistics case as outlined in Section II-C3. 

IV. Numerical Results 

We present numerical results to the schemes discussed in this section. For the numerical results, all the 
noise variances are taken to be equal, i.e. a 2 sd = a 2 si = a 2 d = a 2 . The average SNR of all the links are set 
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Fig. 1: Comparison of SER performance. 



at 0.5, i.e. 7 Qi = 7^ = 0.5, for i. All the shape parameters are taken to be 1 except when indicated. The 
results shown are for a system with 3 relays (m = 3). The relay gain for all relay is taken to be 1. The 
peak individual constraints are set as E r s aax = 3 and E™- ax = 3 for all i. For OPA, E tot is taken to be 5.5 
and for energy-efficiency, j th is taken to be 5 dB. Also, for the energy-efficiency problem, it might be 
the case that due to channel conditions the constraint on the end-to-end SNR cannot be met. In that case 
all transmitting relays and the source transmit at their maximum power. 

Figure 1 shows the comparison of the SER for the different schemes for binary PSK as a function 
of 7 P 8 . Firstly, we compare among the AP and selection schemes for the different cases of CSI at the 
relays. It is evident from Figure 1 that the four OPA AP schemes comfortably provide better performance 
than uniform power allocation (UPA(AP)). Full CSI (AP) gives the greatest gain, as one would expect, of 
around 3.4 dB over UPA (AP) at a SER of 1(T 2 , while partial CSI-a (AP), partial CSI-/3 (AP) and channel 
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statistics (AP) display a gain of around 2.8 dB, 2.6 dB and 2.4 dB, respectively at the same SER. As can 
be noted from the gains of the four OPA schemes over UPA (AP), partial CSI-a (AP) outperforms partial 
CSI-/3 (AP) by a slim margin. However, it must be remembered partial CSI-a (AP) requires knowledge 
of an additional link, namely the source-destination link. Moreover, at a SER of 1CT 2 , all the four OPA 
AP schemes lie within a range of approximately 1.1 dB. Coming to the four OPA schemes for selective 
relaying. As expected, the pattern is similar to the AP case, the four OPA schemes perform better than 
UPA (Sel) 9 . However, the difference in performance for the four OPA selective relaying schemes is large 
as compared to the difference in OPA AP case. Full CSI (Sel) provides a gain of almost 0.8 dB over partial 
CSI-a (Sel) at a SER of 10~ 2 while partial CSI-a (Sel) in turn gives a gain of around 1 dB over partial 
CSI-/3 (Sel). Similarly, Channel Statistics (Sel) trails partial CSI-/3 (Sel) by 1.3 dB. However, Channel 
Statistics (Sel) only provides a gain of 0.2 dB over UPA (Sel). Thus, unlike AP, in selective relaying, 
there is quite a large difference in performance of the four OPA schemes among themselves and a small 
difference in performance among them and their UPA counterparts. Now comparing AP and selective 
relaying, the difference between AP and selective relaying increases with less CSI the system has. In the 
OPA Full CSI case, the difference between AP and selective relaying is only 0.6 dB, however it increases 
with the decrease in channel knowledge and for knowledge of channel statistics only it reaches 2.6 dB. 
Another interesting point to note is that UPA (Sel) outperforms UPA (AP) at low values of j p w . 

This is due to the fact that even though AP has more relays, the total power is the same for both 
AP and selective relaying. For UPA (AP), this power is equally distributed among the 3 relay and the 
source, however, for selective relaying the power is shared between only two nodes and moreover, the 
relay which maximizes the end-to-end SNR is selected. Thus, more power allocated to the relay which 
has better channel conditions and UPA (Sel) performs better than UPA (AP). However, as 7 P increases, 
all the relays see good channel conditions, in general, and the gain of AP is seen. 

The throughputs of all the OPA schemes are shown in Figure 2. All the selective relaying schemes give 
better throughput than all the AP schemes. This is due to the orthogonal distribution of sources in AP. 

'For selective relaying, each case of CSI will have a different UPA scheme because of the different selection criterias. However, as to 
not make the figure cluttered we only show the case for channel statistics. All the other UPA selection schemes show a similar behaviour 
to their corresponding optimal schemes as UPA for channel statistics exhibits with Channel Statistics (Sel). Hence, UPA (Sel) refers to UPA 
for channel statistics. 

I0 It must be kept in mind that UPA (Sel) refers to the case of knowledge of channel statistics only. Hence UPA for selective relaying, 
which has close to Full CSI (Sel performance as mentioned previously), outperforms UPA (AP) at moderate and high values of j p . 
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Fig. 2: Comparison of Throughput. 



Hence, transmitting one packet of information requires m+1 time slots for AP while it requires only 
2 slots for selection. Furthermore, the same pattern is observed for the throughputs as was seen for the 
SER for the CSI case. Partial CSI-a outperforms Partial CSI-/3 and knowledge of only channel statistics 
gives the worst performance. 

Figures 3 and 4 show the performance of AP and selective relaying for the energy-efficiency problem, 
respectively. We will first discuss energy-efficiency for AP relaying by considering Figure 3 and then move 
onto Figure 4 for selective relaying. In the end, we will compare AP and selective relaying. Now, Figure 
3(a) shows the total power consumed, E tot , as a function of 7,, 11 and Figure 3(b) show R 12 as a function 
of 7 S . Figure 3(a) shows that Full CSI (AP) consumes the least amount of energy at low 7 S , however, as 
7s increases, Full CSI (AP) consumes the most amount of energy. This is a little counter intuitive as one 

i2 R = log 2 (l + 7). Hence, R is a linearly increasing function of 7. 
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would expect Full CSI (AP) to give the best performance. However, the analysis is incomplete without 
looking at the achieved end-to-end SNR which is shown by 3(b). The threshold line represents represents 
the value of R for the constraint on the end-to-end SNR, i.e. threshold= log 2 (l + 7 t/l ). As one can see 
from Figure 3(b) that even though Full CSI (AP) consumes the most energy at high 7 S , it also achieves 
the threshold constraint on average. Moreover, at low 7 S , Full CSI (AP) consumes a little less energy 
than all the other CSI cases and achieves the same R on average. The most energy consumed after Full 
CSI (AP) is by Partial CSI-a (AP), however it also provides the second best R after the full CSI case. 
Partial CSI-/3 (AP) and Channel Statistics (AP) achieve the worst R. Selective relaying exhibits a similar 
behaviour in Figure 4 with one major difference. At high Partial CSI-a (Sel) achieves a greater R 
than Full CSI (Sel) which approximately achieves the threshold value. However, this comes at increased 
energy consumption as shown by Figure 4(b) due to the mismatch in instantaneous values and the average 
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Fig. 5: Effect of shaping parameter. 



value. 

Now comparing AP relaying and selective relaying, all the selective relaying cases consume less energy 
than their AP counterparts for low values of 7 S which is understandable as when the channel conditions 
are bad, AP will consume 12 energy units while selection only 6 due to only one relay being active. 
Consequently, the selective relaying CSI cases achieve a lower R than their corresponding AP cases. 
However, this pattern continues even for at high 7 S where AP would expect to consume less energy. The 
reason for this is again the case when the channel conditions are bad and the constraint on the end-to- 
end SNR cannot be met. Even though bad channel conditions occur rarely for at high 7 S , but, in these 
instances, AP consumes twice the power of selective relaying which pushes the average power consumed 
for AP a little above than that of selective relaying. Therefore, it seems that for energy-efficiency, selective 
relaying is a better option than AP even at high 7 S . 



22 

The effect of the shaping parameter A; 13 is shown in Figure 5 14 for j p =0 dB. As can be seen from 
Figure 5, a higher shape parameter means better performance for all the CSI cases. However, this does 
not effect the difference in the three CSI cases. 

V. Conclusions 

We have studied power allocation to maximize the end-to-end SNR under a total power constraint and to 
minimize the the total power consumed while maintaining the end-to-end SNR above a required threshold 
for a fixed-gain AF relay network. We have studied both problems for the relay network operating in AP 
mode where all the relays participate in signal forwarding and operating in selection mode in which only 
the selected relay forwards the signal to the destination. Furthermore, we have also considered the cases 
of full CSI, partial CSI and knowledge of only channel statistics at the relay for both modes of operation 
and for both optimization problems. We demonstrate the gain achieved by allocating power optimally over 
UPA. We also give insight into the performance of the system for both problems, for both AP relaying and 
selective relaying and for all three cases of CSI at the relays. Additionally, we also develop inequalities in 
Appendix B which may prove to be useful in future works. For these reasons, we believe that our work 
is a valuable contribution to the already available literature on power allocation strategies for fixed-gain 
AF relays. 

Appendix A 

Writing down the objective function 

The objective function, in general, not convex and concave. However, as we show below, it is concave 
(its negative is convex) for the domain we are interested in. 
Define vector E as 

E=[E 8 E 1 E 2 ■■■ E m } T E y 0. (43) 

l3 k ai = kp t — k for all i 

14 We only show the results for AP as selection also displays similar results and behaviour and they are omitted so as to not clutter the 
plot and make it difficult to follow. 
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Now let us define 



/(E)=[100...0]E = E S g(E)=J2^-J2 



Both / and g are positive and increasing on their domain. For / to be concave 



(44) 



f(0x + (l-0)y)>0f(x) + (l-0)f(y), 



(45) 



where < 9 < 1. The left hand side (LHS) in the above is 9x 1 + (1 — 9)y 1 and the right hand side (RHS) 
is equal to 9x\ + (1 — 9)y\. As the LHS is equal to the RHS, / is concave. To show that g is concave, 
forming the Hessian 





(aiBift+Cl) 3 









(a 2 £ 2 /3 2 +C2) 3 

















(tlmJ?m(Sm+Cm) :! , 



(46) 



As the eigenvalues of H s are non-negative, H g is negative semi-definite, and hence g is concave. Now let 
us define 



fc(E ) = / ( EME) = £ ,( Eai -g-^L_ f _ 



(47) 



,i=0 



For /i to be concave (— h to be convex) 



h{9x + (1 - %) > 9h(x) + (1 - 6>)/i(y). 



(48) 



Therefore for concavity we have to show 



A < 0, 



(49) 



where 



A = 9(fg)(x) + (1 - 9){fg){y) - (fg)(9x + (1 - 9)y) 



(50) 
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As / and g are both positive and concave 

(fg)(9x + (1 - 9)y) > (9f(x) + (1 - 9)f(y))(9g(x) + (1 - 9)g{y)) (51) 

Substituting (51) in the expression of A one has 

A < 9{fg){x) + (1 - 9){fg){y) - (9f(x) + (1 - 9)f{y)){9g{x) + (1 - 9)g{y)). (52) 
A < 9f{x)g{x) + (1 - 9)f{y)g{y) - 9 2 f{x)g{x) - (1 - #) 2 /(2/M?/) 

(53) 

-e(l-0)f(x)g(y)-e(l-e)f(y)g(x). 

After some manipulation 

A<0(l-0)D(z,2/), (54) 

where D(x,y) = (f(x) — f(y))(g(x) — g{y)). If D(x,y) < 0, then the proof of concavity is complete. 
For optimal power allocation, if f(x) > f(y) then ^2 iey Ei > ^2 iex Ei from the total power constraint. 
Therefore, with power allocation, g(y) > g(x), as it g(.) is a concave functions of the relay powers, 
implying D(x,y) < and hence, concavity. Similarly, if g(x) > g(y), then power allocation means that 
J2iex Ei > J2ie y Ei, which in turns means f(y) > f(x). Thus, D(x, y) < and the objective function is 
concave. 

Appendix B 

Writing down the objective function 

It is obvious that it is a concave function of £" s . To check concavity with respect to Ej, we find the second 
derivative 

= - k . it > k ^ (v,^), (55) 
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Fig. 6: Plot of S(v,k 0j ). 

where v = — J^- and 

S(v, k Pj ) = -(2 + k Pj )v~ k ^ +1 - v~ k ^ +2 + T (1 - k Pj ,v) e v (v 2 + 2k Pj v + 2v + k\ + (56) 

Therefore, for concavity with respect to £j, we have to show 

> v>0, kfr >0. (57) 

Figure 6 shows the plot of £>(.). It is apparent from Figure 6 that S'(.) > 0. However, we did not show S(.) 
for very low and high values of v and k^). The is due to the fact that for low values fo the parameters, 
S(.) has a very high value exceeding 10 20 as can be seen from the pattern of Figure 6. At very high 
values of v and k^), the values of S(.) become very low even falling below 1CT 18 after which software 
packages like MATLAB can't accurately characterize the number on plots. However, the value does not 
become negative at any point. Therefore, the objective function is concave with respect to Ej. Now using 
a similar argument as in Appendix A, it can be shown that the objective function jointly concave with 
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respect to E s and Ej. 

For the special case of Rayleigh fading, the second derivative becomes 



= ^ 7aj I -3 - + ( -^—] e a ^i ^ _ 

dE) S ajE*% \ ' ajErfp. \ajEj* J "V'j'X 




(58) 
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In this case, the proof of concavity establishes the relationship 

E 1 (v) e v (v 2 + 2 + Av) - v - 3 > v > 0. (59) 
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